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Abstract. We have investigated magnetically driven superlu- 
minal jets originating from rotating black holes. The station- 
ary, general relativistic, magnetohydrodynamic wind equation 
along collimating magnetic flux surfaces has been solved nu- 
merically. Our jet solutions are calculated on a global scale 
of a spatial range from several to several 1000 gravitational 
radii. Different magnetic field geometries were investigated, 
parameterized by the shape of the magnetic flux surface and 
the magnetic flux distribution. For a given magnetic flux sur- 
face we obtain the complete set of physical parameters for the 
jet flow. In particular, we apply our results to the Galactic su- 
perluminal sources GRS 1915+105 and GRO 1655-40. Moti- 
vated by the huge size indicated for the Galactic superluminal 
knots of about 10 9 Schwarzschild radii, we point out the pos- 
sibility that the jet collimation process in these sources may be 
less efficient and therefore intrinsically different to the AGN. 
, Our results show that the observed speed of more than 0.9 c 
fSJ ' can be achieved in general by magnetohydrodynamic acceler- 
ation. The velocity distribution along the magnetic field has a 
r ~ l ' saturating profile. The asymptotic jet velocity depends either on 
, the plasma magnetization (for a fixed field structure) or on the 

f ' magnetic flux distribution (for fixed magnetization). The dis- 

(f) , tance where the asymptotic velocity is reached, is below the ob- 
t— I ' servational resolution for GRS 1915+105 by several orders of 
<0 magnitude. Further, we find that highly relativistic speeds can 
be reached also for jets not emerging from a region close to the 
black hole, if the flow magnetization is sufficiently large. The 
i-^ ' plasma temperature rapidly decreases from about 10 10 K at the 
^p" 1 , foot point of the jet to about 10 6 K at a distance of 5000 gravita- 
O ' tional radii from the source. Temperature and the mass density 
+3 , follow a power law distribution with the radius. The jet mag- 
^ netic field is dominated by the toroidal component, whereas the 
• • , velocity field is dominated by the poloidal component. 
> ■ 

• i-H 

, Key words: Accretion, accretion disks Black hole physics - 
;_j ' MHD - Stars: mass loss - ISM: jets and outflows - Galaxies: 
\ jets 
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1. Introduction 

1.1. Relativistic jets and Galactic superluminal motion 

Apparent superluminal jet motion originating in the close envi- 
ronment of a rotating black hole is observationally indicated for 
two classes of sources concerning mass and energy output. One 
class is the family of radio loud active galactic nuclei (hereafter 
AGN). In the AGN standard model highly relativistic jet mo- 
tion is explained by magnetohydrodynamic processes in a black 
hole - accretion disk environment (for a review see Blandford 
1990). Jets are magnetically accelerated and possibly also col- 
limated by magnetic forces. However, the detailed interaction 
process of the magnetized black hole - accretion disk system 
which is believed to lead to the ejection of high velocity blobs 
is not yet fully understood. 

The other class are galactic binary systems for which 
radio observations have also detected superluminal mo- 
tion (see reviews of Fender 2000 or Greiner 2000). The 
two most prominent examples are the high energy sources 
GRS 1915+105 (Mirabel & Rodriguez 1994) and GRO 1655- 
40 (Hjellming & Rupen 1995; Tingay et al. 1995). The de- 
projected jet speed of both sources is <; 0.9 c and surpris- 
ingly similar, although for GRS 1915+105 also a higher veloc- 
ity component has been observed recently (Fender et al. 1999). 
GRO 1655-40 is a binary consisting of a 7.02 ±0.22 M black 
hole and a 2.3 Mq F-subgiant (Orosz & Bailyn 1997) at a dis- 
tance of 3 kpc. GRS 1915+105 is at 10-12 kpc distance (Fender 
et al. 1999), but the component masses of the presumed binary 
are not known. Order of magnitude estimates based on X-ray 
variability and QPO properties range from 10-80 Mq (Morgan 
etal. 1997, Greiner etal. 1998). As for the AGN jet sources, ob- 
servational evidence for a black hole - accretion disk system is 
found also for the Galactic superluminal sources. Observations 
have also indicated that accretion disk instabilities may be re- 
lated to jet ejection (Greiner et al. 1996, Belloni et al. 1997, 
Mirabel etal. 1998). Therefore, the jet formation process for 
extragalactic jets and their Galactic counterparts may be the 
same, although the mechanism that accelerates and collimates 
the GRS 1915+105 ejecta is yet unclear (Rodriguez & Mirabel 
1999). 
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Optical polarization measurements have been obtained for 
the microquasar GRO J1655-40 (Scaltriti et al. 1997, Gliozzi 
et al. 1998). The polarization angle is approximately parallel to 
the accretion disk plane. The amount of polarization has been 
found to vary smoothly with the orbital phase, being smallest 
at binary phase 0.7-0.8. It has been noted that the occasionally 
observed X-ray dips occur at the same phase interval (Ueda et 
al. 1998, Kuulkers et al. 1998) suggesting that it may be related 
to either a thickening of the disk rim at the impact site of the 
accretion stream from the companion or the overflow of this 
stream above/below the disk. The orbital polarization modula- 
tion rules out a synchrotron origin in the jet, and implies the 
presence of electron scattering plasma above the accretion disk 
which is asymmetrically distributed or asymmetrically illumi- 
nated. The existence of such scattering plasma is consistent 
with the interpretation of the iron features as observed with 
ASCA as absorption lines and edges in a thick, cool torus of 
column N H > 10 23 cirr 2 (Ueda etal. 1998). 

The relativistic speed observed for the Galactic superlumi- 
nal sources (~ 0.9 — 0.98 c de-projected) corresponds to a bulk 
Lorentz factor of 7 = 2 — 5 although this number is not very 
accurate (e.g. Fender et al. 1999). Therefore, for any theoretical 
investigation of these objects at least special relativity has to be 
taken into account. If the superluminal motion originates close 
to a black hole, also general relativistic effects may become 
important. 

The ejection of matter itself is not a stationary process. 
In GRS 1915+105 also repeated emission of knots is observed 
(Rodriguez & Mirabel 1999). X-ray and radio observations 
suggest that a wide range of ejected mass and ejection fre- 
quency is possible. 

Though the galactic jet sources are nearby, they are not bet- 
ter resolved spatially because the distance ratio between AGN 
and microquasars is smaller than their mass ratios. Neverthe- 
less, an important implication may also come from the ob- 
served size of the superluminal knots which are observation- 
ally resolved. In the case of GRS 1915+105 the characteristic 
dimension of the 'jet' is 35 mas, equivalent to 7 x 10 15 cm at a 
distance of 12.5 kpc (Rodriguez & Mirabel 1999). We empha- 
size that such a knot size corresponds to ~ 10 9 Schwarzschild 
radii for R s = 1.5 x 10 6 (M /5 M Q )cm! This is a huge factor 
and may be in distinctive difference to the AGN jets. Similarly, 
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the VLBA data show the core as a collimated jet down to a dis- 
tance of 10 AU from the central source with an opening angle 
of < 10° (see Mirabel & Rodriguez 1999) corresponding to 
10 7 (M/5Mq ) Schwarzschild radii. The length of the radio jet 
is about 100 AU. 

However, when interpreting the observed emission region, 
one has to keep in mind that this region may not represent the 
jet flow itself, but some part of another, larger, structure. For 
example, in some extragalactic jet sources there is indication 
that the knots travel along helical trajectories, believed to be 
prescribed by a large-scale helical magnetic field of an almost 
cylindrically collimated jet (Zensus et al. 1995; Camenzind & 
Krockenberger 1992). 

In GRO 1655-40 the motion of the radio knots is com- 
plicated and requires (at least) precession between different 
ejections (Hjellming & Rupen 1995). The knot structures in 
GRS 1915+105 remained fixed implying that the whole knot 
moves with the same speed without spatial diffusion and with 
an axial velocity profile more or less constant. 

Based on minimum energy arguments and only relativis- 
tic electrons responsible for the synchrotron radiation in the 
knots of GRS 1915+105 , Rodriguez & Mirabel (1999) derive 
a magnetic field strength of about 50 mG to 7 mG, the decrease 
resulting from the expansion of the knot. They also estimate 
the rest mass of a knot of > 10 23 g, and together with (steady) 
photon luminosity of ~ 3 x 10 38 ergs -1 , exclude radiation as 
driving mechanism for the knots. 

1.2. Theory of magnetic jets 

From the introductory remarks it is clear that a quantitative 
analysis of superluminal motion must take into account both 
magnetohydrodynamics (hereafter MHD) and (general) rela- 
tivity. 

The first theoretical formulation of the electromagnetic 
force-equilibrium in Kerr space-time around rotating black 
holes was given by Blandford & Znajek (1977) and Znajek 
(1977), who discovered the possibility of extracting rotational 
energy and angular momentum from the black hole electro- 
magnetically. 

Camenzind (1986, 1987) formulated a fully relativistic sta- 
tionary description of MHD flows, basically applicable to any 
field geometry. The structure of such collimating jet magneto- 
spheres in the case of Kerr space time was presented by Fendt 
(1997). Solutions of the so-called wind equation in Kerr ge- 
ometry (see below) considering the stationary plasma motion 
along the magnetic field were obtained by Takahashi et al. 
(1990), however, mainly discussing the accretion flow onto the 
black hole. 

While the asymptotic structure of the propagating jets be- 
comes more and more understood with the help of time- 
dependent magnetohydrodynamical, also relativistic, simula- 
tions (e.g. Nishikawa et al. 1997; Mioduszewski et al. 1997; 
Hardee et al. 1998), the process of jet formation itself and the 
collimation of the outflow region is a task still too complex for 
numerical simulations. The involved length scales and gradi- 



ents require a high resolution in grid size and time stepping. 
Koide et al. (1998) were first to perform general relativistic 
MHD simulations of jet formation close to the black hole. In 
their model, the interaction of an initially cylindrical magnetic 
field with a Keplerian accretion disk results first in an inflow 
of matter towards the black hole. This accretion stream inter- 
acts with the hydrostatic corona around the black hole giving 
rise to a relativistic gas pressure driven jet. At larger radii a 
magnetically driven wind is initiated from the accretion disk. 
The simulations were performed for less than two rotations of 
the inner disk (corresponding to less than 0.02 rotations of the 
disk at the outer edge of the grid). Although these results of 
the first fully general relativistic MHD simulations look indeed 
very exciting, some objections can be raised about the under- 
lying model. The initial condition applied is that of a hydro- 
static corona around a black hole, an assumption which is not 
compatible with the boundary of a black hole horizon. Such a 
configuration is not stable and will immediately collapse. Re- 
cently, the authors extended their work applying an initial coro- 
nal structure in steady infall surrounding a non-rotating black 
hole (Koide et al. 1999). They find a two-layered jet consisting 
of a magnetically driven jet around a gas-pressure driven jet. 
In addition, Koide et al. (2000) considered the quasi-steady in- 
fall of the corona around a Kerr black hole. They find that jet 
formation seems to differ for co-rotating and counter-rotating 
disks. The jet ejection tends to be easier in the latter case with 
a jet origin much closer to the hole. Also, a new feature of an- 
other magnetically driven (though sub-relativistic) jet appears 
within the gas-pressure driven jet. The computations were last- 
ing over a few inner disk orbits. Therefore, the observed events 
of mass ejection could still be a relict of the initial condition 
and may not be present in the long-term evolution. Clearly, it 
would be interesting to perform the Koide et al. simulations 
for a longer time and look whether the mass ejection continues 
over many disk orbits, whether the simulation evolves into a 
final stationary state (as e.g. in Ouyed & Pudritz 1997, Fendt 
& Elstner 2000), or whether the jet formation retains its un- 
steady behavior which could explain the emission of superlu- 
minal knots observed in the relativistic jets. 

1.3. Aim of the present study 

In this paper, a stationary magnetic jet flow along a given mag- 
netic flux surface is investigated in the context of general rela- 
tivity. Due to the stationary approach, we cannot treat any time- 
dependent phenomena. Our emphasis is to trace the large scale 
behavior of the flow from it's origin close to the black hole 
to large distances. This is an essential point in particular for 
the Galactic superluminal sources because of the possible huge 
spatial extension of the jets compared to the central black hole. 
The stationary model allows for a global treatment of the jet 
flow, i.e. an investigation over a large range of magnitudes for 
density and magnetic field strength. This is not yet feasible with 
time-dependent MHD codes presently available. In particular, 
we address the following topics. 
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- For a given geometry of the magnetic field, which are the 
resulting jet dynamical parameters as velocity, density or 
temperature? 

- How important are the effects of general relativity? Does 
the superluminal flow indeed originate very close to a black 
hole? 

- From the investigation of different field geometries we ex- 
pect some hints to the jet opening angle and the length scale 
of the collimation process. 

The structure of this paper is as follows. In Sect. 2, basic equa- 
tions for relativistic magnetospheres are reviewed in the context 
of Kerr metrics. In Sect. 3, the model underlying our numerical 
calculations is discussed. We present our numerical results in 
Sect. 4 and discuss solutions with different geometry and jet 
parameters. We summarize our paper in Sect. 5. 

2. Description of a MHD flow in Kerr metric 

Under the assumptions of axisymmetry, stationarity and infi- 
nite conductivity, the MHD equations reduce to a set of two 
basic equations describing the local force-balance across the 
field and along the field (for references, see, e.g., Blandford 
& Znajek 1977; Thorne et al. 1986; Camenzind 1986, 1987; 
Okamoto 1992; Beskin & Pariev 1993, Beskin 1997). 

The trans-field or Grad-Shafranov equation determines the 
field structure, whereas the wind equation describes the flow 
dynamics along the field. Due to the stationarity assumption, 
certain conservation laws apply. The total energy density, the 
total angular momentum density, the mass flow rate per flux 
surface and the iso-rotation parameter are conserved quanti- 
ties along the surfaces of constant magnetic flux (Camenzind 
1986). 

In this paper the motion of a magnetized plasma is calcu- 
lated from the wind equation. The plasma moves along a pre- 
scribed axisymmetric magnetic flux surface which originates in 
a region close to a rotating black hole. 

2.1. Space-time around rotating black holes 

The space-time around a rotating black hole with a mass M and 
angular momentum per unit mass a is described using Boyer- 
Lindquist coordinates with the line element 
ds 2 = a 2 dt 2 - d, 2 [d4> - udt) 2 - {p 2 /A) dr 2 - p 2 d6 2 , (1) 
where t denotes the global time, cf> the angle around the axis of 
symmetry, r, 6 similar to there flat space counterpart spherical 
coordinates, and where geometrical units c = G = 1 have been 
applied (see Appendix A for further definitions). The horizon 
of the Kerr black hole is located at rjj = M + V M 2 — a 2 . 
We will normalize all radii to gravitational radii r g = rn(a = 
M) = Al. The angular velocity of an observer moving with 
zero angular momentum (ZAMO) is u> — (d(j>/dt)zAMO, cor- 
responding to the angular velocity of the differentially rotating 
space. The lapse function is a = (dr/dt)zAMO describing the 
lapse of the proper time t in the ZAMO system to the global 
time t. 



2.2. Description of the electromagnetic field 

In the 3+1 split of Kerr space time (Thorne et al. 1986) the elec- 
tromagnetic field B, E, the current density j, and the electric 
charge density p c can be described very similar to the usual ex- 
pressions, if measured by the ZAMO's according to the locally 
flat Minkowski space. These local experiments then have to be 
put together by a global observer for a certain global time us- 
ing the lapse and shift function for the transformation from the 
local to the global frame. 

With the assumption of axisymmetry a magnetic flux sur- 
face can be defined measuring the magnetic flux through a loop 
of the Killing vector m = a) 2 V(/>, 



(2) 



If 1 

i$>(r, 0) = — -Bp ■ dA , B p = — V* A m, 

Z7T J UT 

corresponding to the magnetic flux through an area ir(r sin 6) 2 
around the symmetry axis (in the limit of Minkowski space). 

With the assumption of a degenerated magnetosphere, 
| \B\ 2 - \E\ 2 \ » \E ■ B\ ~ an 'angular velocity of field 
lines' can be derived from the derivative of the time component 
of the vector potential tt F = £1 F (*) = -2itc(dA / d^). We 
will denote this quantity with the term 'iso-rotation parameter' . 



2.3. The wind equation 

It has been shown that a stationary, polytropic, general relativis- 
tic MHD flow along an axisymmetric flux surface >J/(r, 6) can 
be described by the following wind equation for the poloidal 
velocity u p = "fv p /c, 



1 



k k 2 + a m 2k 2 M 2 - k 4 M 4 



(k 



M 2 y 



(3) 



where 

h = 



.933^1 + 2 303^F 

L 



5oo, 



k 2 = 1 - F — , 
E 



( „ L L 2 \ . . 2 , 

«4 = ~ I .933 + 2.903"^ + 500-^2" I / (.903 ~ .900533) 

(Camenzind 1986, Takahashi et al. 1990). The Alfven Mach 
number M is defined as M 2 = 4/j,nu 2 /B 2 , with the proper 
particle density n, the specific enthalpy /i, and a poloidal mag- 
netic field i? p = -Bp/(g o)+go30F), rescaled for mathematical 
convenience. The quantity a m stands for the sign of the metric 
(we have chosen a m = — 1, see appendix A). For a polytropic 
gas law with the index T = n/m, the wind equation (3) can be 
converted into a polynomial equation, 

2n+2m 

A ?; (x;*,$;ft F ;£,L, ( 7*)wp /m = 0, (4) 

i=0 

(Camenzind 1987, Englmaier 1993, Jensen 1997), where the 
coefficients A4 are now defined as functions of the normal- 
ized cylindrical radius x = R/r g (see Appendix B). The shape 
of the axisymmetric magnetic flux surface is prescribed as 



function z(x; i b). The flux function $ 



gB p describes the 



opening of the flux tube. The faster $ decreases the faster mag- 
netic energy is converted into kinetic energy. We define the di- 
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mensionless magnetization parameter [] at the 'injection' point 
x+ following Takahashi et al. (1990), 

a * = A V' (5) 

measuring the Poynting flux in terms of particle flux I p = 
\J—g nu p , where m p is the particle mass (here the proton 
mass). The magnetization determines the maximum energy 
available for plasma acceleration and thus determines also the 
asymptotic poloidal velocity. The other wind parameters are to- 
tal energy density E, total angular momentum L, and the iso- 
rotation parameter Sip- The non relativistic limit of Eq. (4) has 
been solved numerically by Kudoh & Shibata (1995, 1997). 

We choose the polytropic index T = 5/3 for a hot rela- 
tivistic proton-electron plasma (a hot electron-positron plasma 
would imply T = 4/3). Then, at each radius x the polynomial 
equation (4) has 2n + 2m = 16 solutions. Some of these math- 
ematical solutions have no physical meaning, e.g. because u p 
is negative. The remaining physical solutions form a bunch of 
different curves in the u p (a;)-diagram representing different so- 
lution branches (see our solution SI in Appendix C, Fig. CI). 
The unique branch of the 'wind' solution starts at a small radius 
with small velocity continuing outwards with increasing ve- 
locity. For an other parameter choice also 'accretion' branches 
can be found, starting from a large radius with small velocity 
and continuing inwards with increasing velocity (not shown in 
Fig. CI). 

However, not for all parameters E, L, a there exist physi- 
cal solutions which are continuous functions of x and there- 
fore defined along the whole flux surface. It is well known that 
at the magnetosonic points the wind equation (3) becomes sin- 
gular (see Camenzind 1986, Takahashi et al. 1990). Regularity 
of the solution requires a flow velocity equal to the speed of 
the MHD waves in order to obtain a smooth (self-consistent) 
transition at the magnetosonic points. In order to match astro- 
physical boundary conditions we fix the following parameters, 

- the 'injection' radius, a;*, the location where the matter 
couples to the magnetic field. This radius also determines 
the iso-rotation parameter ftp, 

- the 'injection' velocity u p * = u p (x i ,), defining the initial 
kinetic energy. 

- the Alfven radius x&, which fixes the total angular momen- 
tum of the flow. 

The critical wind solution for a given flux surface can then be 
found by varying the flow parameters in Eq. (4). Due to numer- 
ical convenience, we vary 

- the sound speed c s * at the injection radius, defining the ini- 
tial density (or gas pressure and temperature), 

1 Note that this definition for the magnetization varies from the 
original Michel magnetization parameter <tm = &m/A-k fmcR^, 
where <E>m is the magnetic flux, /m the mass flux and _Rl the light 
cylinder. Usually, the general relativistic equations are normalized to 
the gravitational radius, whereas the special relativistic equations are 
normalized to the light cylinder 



Fig. 1. Model geometry applied for our numerical calculations. 
The poloidal field structure is prescribed as magnetic flux sur- 
faces with different opening angle. The flux surfaces have dif- 
ferent foot point radii along the accretion disk (not visible). The 
central source is a black hole implying that general relativistic 
effects have to be taken into account. The toroidal field follows 
from the solution of the wind equation. 

- the magnetization parameter at the injection point er* = 

$2/(47rm p / p ,). 

In turn, the condition of a regular flow at the magnetosonic 
points fixes the sound speed and magnetization and, thus, jet 
mass flow rate and temperature. 

3. The model assumptions 

3.1. The model in general 

Observationally the jet phenomenon of AGN, young stellar ob- 
jects and microquasars is always connected to the signatures of 
an accretion disk. We therefore assume a similar disk-jet sce- 
nario for the jet formation in Galactic superluminal jet sources. 
In general our model geometry follows the standard model for 
jet formation in AGN (cf. Blandford 1990). 

Two typical length scales enter the problem, (i) The gravita- 
tional radius r g measures the influence of gravity on the metric, 
(ii) The asymptotic light cylinder describes the influence of 
rotation on the electrodynamics. 

3.2. The central black hole 

The black hole mass and angular momentum determine the ge- 
ometry of space. Since we use dimensionless equations nor- 
malized to the gravitational radius, our results scale with the 
mass of the black hole. For parameter estimates we assume 
a black hole mass of 5 Mq which is about the value inferred 
for the galactic superluminal sources. The angular momentum 
a as the other black hole parameter is not known for any of 
the relativistic jet sources. Interpretation of the high effective 
temperatures of the accretion disk as well as the stable QPO 
frequency (as Thirring-Lense effect) suggests that a ^20.9 for 
GRS 1915+105 and GRO 1655-40 (Zhang etal. 1997). Theo- 
retically, one may expect a rapidly rotating black hole because 
of angular momentum conservation during the collapse and 
also accretion of angular momentum from the accretion disk 
(King & Kolb 2000). Here, we choose a = 0.8, a value which 
is not extreme, but clearly different to Schwarzschild metric. 
The rotation rate of the black hole is defined as fin = w(^h) = 
a/ (2Mr-n). The Kerr parameter a does not influence the solu- 
tion of the wind equation directly. However, for rotating black 
holes the marginally stable orbit r ms is closer to the horizon, 
r ms = 6r g for a = and r ms ~ r g for a ~ 1 (This is the case 
for a co-rotating disk. For a retrograde disk rotation r ms ~ 9r g 
for a ~ 1). Therefore, assuming that the jet magnetic field is 
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anchored just at the marginally stable orbit, for a rapidly ro- 
tating black hole the maximum angular velocity of the jet foot 
points increases by a factor of 6 3 / 2 /2 = 7.4. Correspondingly, 
the light cylinder radius of the jet moves inward by the same 
factor. 

In addition to the well-known special relativistic light cylin- 
der, the differential rotation of the space u> leads to the forma- 
tion of a second light surface. At this position the 'rotational ve- 
locity' of the field lines relative to the ZAMO equals the speed 
of light (see Blandford & Znajek 1977). The position of the 
two light surfaces <1>l is defined by u) 2 = (±qc/(Of — w)) 2 , 
where the + (— ) sign holds for the outer (inner) light surface 
with fip > uj (f2p < oj). However, these light surfaces have 
no direct implication for the MHD flow. In the limit of a strong 
magnetization, the MHD Alfven surfaces (for inflow and out- 
flow) approach the corresponding light surfaces. 

3.3. The accretion disk 

X-ray observations of GRS 1915+105 detected strong inten- 
sity variations indicating major instabilities of an accretion disk 
(Greiner et al. 1996). Belloni et al. (1997) find that the highly 
variable X-ray spectra could be explained if the inner disk is al- 
ternatively removed and replenished due to a thermal-viscous 
instability. Simultaneous X-ray and infrared observations of 
GRS 1915+105 revealed evidence for a disk-jet interrelation 
(Eikenberry et al. 1998, Mirabel etal. 1998). The observed 
flares in the X-ray and IR bands have a consistent offset de- 
lay of ~ 30 min indicating an origin from the same event. 

The accretion rate in GRS 1915+105 and GRO J1655-40 
can be determined from the observed X-ray luminosities (e.g. 
Greiner etal. 1998). Depending on the chosen efficiency (5% 
in non-rotating versus 42% in maximally rotating black holes) 
the accretion rate ranges between 1 — 9 x 10 _7 MQyr _1 (GRS 
1915+105) and 0.8 - 7 x lO~ 8 M yr- 1 (GRO J1655-40), re- 
spectively. 

From the theoretical point of view an accretion disk sur- 
rounding the black hole is the essential component concerning 
magnetic jet formation. It is considered to be responsible for 
the following necessary ingredients for jet formation, propaga- 
tion, and collimation. 

- The generation of the magnetic field. In contrast to stel- 
lar jets the magnetic field of jets from black holes cannot 
be supplied by the central object but has to be generated 
by the surrounding accretion disk. Dynamo action in gen- 
eral relativistic accretion disks were discussed by Khanna 
& Camenzind (1996a, 1996b) and Brandenburg (1996). 

- The mass loading of the jet. The accreting material be- 
comes partly diverted into the jet. Evidently, no mass out- 
flow is possible from the black hole itself, in difference to 
a stellar wind. The (non-relativistic) self-similar accretion- 
ejection mechanism was investigated by Ferreira (1997). 

- The electric current system. Differential rotation of the disk 
is also responsible for driving the poloidal electric current 
system in the jet magnetosphere. Such a current extracts 



angular momentum from the disk and eventually allows for 
mass accretion into the central object. 

3.4. Model parameters for the wind motion 
3.4.1. The magnetization parameter 

The leading parameter for the wind solution along a 
fixed poloidal field is the magnetization parameter (5). Re- 
normalization to astrophysical units gives 

RA 4 

— (6) 



47rm p / p * cMjetWr 2 cM jct (V)\r g 

where Afj Ct (4') ~ 47rra p n*cu p *i? 2 is the jet mass flux en- 
closed by an area of radius R+. A first order estimate of 
the magnetization can be derived from the disk equipartition 
field strength. Then, with a reasonable assumption on the jet 
mass flow rate related to the disk accretion rate, this gives the 
jet magnetization. Although the equipartition field strength is 
model-dependent, the different models (e.g. either advection 
dominated disk or standard disk, either Kramer's opacity or 
Thomson scattering) give rather similar results. A self-similar 
advection dominated disk model with the accretion rate M acc 
gives 



/ M 



B cq ~ 2.5 x 10 y G a vi i 



(7) 



where M E = 1.1 x 1O" 7 (M/5M ) Mgyr" 1 is the Edding- 
ton luminosity and a v i s is the viscosity parameter (see e.g. 
Narayan et al. 1998). In comparison, an optically thin standard 
accretion disk with Thomson opacity gives B cq ~ V8^P = 
y/8naT*/3 ~ 1.8 X 10 8 Ga; ib 1/2 (A//5M Q )- 1 /2( i? / rg )-3/4 
(see Blandford 1990). Note that these estimates are only valid 
within certain limits of the accretion rate and the disk radius. 
Using the advection dominated disk model equipartition field 
strength, we obtain the following estimate for the magnetiza- 
tion at the injection radius, 



1 



M 



1 jet 



3/2 



(8) 



<r*(#) = 16 . 

a vis V 5M J y M acc / v ' g 

A comparison with the original Michel magnetization parame- 
ter % must take into account a factor (r g /R^) 2 . The magne- 
tization parameter derived from the field distribution in a stan- 
dard accretion disk model (see above) will give a similar result. 
We emphasize that we do not 'apply' a certain disk model (e.g. 
the ADAF model) in our computations. However, a comparison 
in the context of accretion disk theory just puts our wind pa- 
rameters on a safer ground. Note, that neither the ADAF model 
nor the standard disk model takes into account the influence 
of magnetic fields. Moreover, the ADAF estimates as cited in 
Eq. (7) rely on the self-similar assumption. Compared to the 
standard disk, by definition, the matter in the ADAF disk would 
be rapidly advected possibly influencing also the wind ejec- 
tion. However, such a detailed treatment is beyond the scope 
of this paper and may only be considered in numerical simu- 
lations investigating the disk-jet interaction itself (Koide et al. 
1998, 1999, 2000) 
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Table 1. Comparison of leading parameters for the wind solution. Magnetic flux distribution shape of the poloidal field 

line z(x), iso-rotation parameter Op, sound speed at the injection radius c s *, magnetization at the injection radius er*, cylindrical 
Alfven radius xa, cylindrical injection radius x*, total energy E, normalized to m p c 2 normalized total angular momentum 
L = L/E, asymptotic velocity u poo , and angular momentum parameter of the black hole a. Other Parameters are: T = 5/3, 
up* = 0.006 (S3-S9), u pi < = 0.17 (S3q, S3u2), u p + = 0.21 (S3u3) 



prescribed calculated 





$/$ 


t 




z{x) 




Of 




at 


XA 


,r. 


E 


L 


Upoo 


a 


S3 


~ 1 




0.1 


[x 


- x 


6/5 


0.035 


0.05165 


979.4 


22.931 


8.3 


2.7887 


20.04 


2.531 


0.8 


S3c2 


~ 1 




0.1 


X 


- x 


3/2 


0.035 


0.0529 


1356 


22.931 


8.3 


2.764 


19.95 


2.58 


0.8 


S4 


~ x~~ 


1/2 


0.1 


[x 


- x 


6/5 


0.035 


0.049 


2380 


22.931 


8.3 


2.7879 


20.04 


2.60 


0.8 


S4b 


~ x~ 


1/2 


0.1 


[x 


- x 


6/5 


0.014 


0.0390 


14680 


57.0 


15.3 


2.6730 


47.07 


2.48 


0.8 


S9 


~ x~ 


1/2 


0.1 


[x 


- x 


3/2 


0.035 


0.05165 


2777 


22.92 


8.3 


2.7572 


19.93 


2.57 


0.8 


S3q 


~ 1 




0.1 


[x 


- x 


6/5 


0.14 


0.31 


480 


5.83 


3.3 


8.917 


6.616 


8.48 


0.8 


S3u 


~ 1 




0.1 


[x 


- x 


6/5 


0.14 


0.27 


100 


5.33 


3.3 


3.16 


5.69 


2.96 


0.8 


S3u2 


~ 1 




0.1 


[x 


- x 


6/5 


0.14 


0.27 


82.5 


5.33 


3.3 


4.66 


6.35 


4.55 


10~ 8 


S3u3 


~ 1 




0.1 


[x 


- x 


6/5 


0.14 


0.27 


205.7 


5.33 


3.3 


4.65 


6.35 


4.48 


10" 8 



Fig. 2. Projected magnetic flux surface. Shape of the poloidal 
field line / flux surface as function z(x) for the solutions S4 
(and S4b, S3, SI) and S9 (and S3c2). 



3.4.2. The magnetic field distribution 

The normalized magnetic field distribution is prescribed by 



the shape of the field line, z(x), 



— the magnetic flux distribution, <&(x) = <&(x) yj— g/(p 2 A). 

We apply different functions for z(x) and $(a;) in order to in- 
vestigate the influence of collimation, rotation and magnetic 
flux distribution on the acceleration of matter. One example is 
z(x) = 0.1(x — xo) 6 ^ 5 describing an almost conical surface 
with only a slight collimation (see Fig. 2). Here, xq defines the 
intersection of the field line with the equatorial plane, with xq 
somewhat smaller than ir*. The idea behind this choice is that 
the matter is expected to couple to the jet magnetic field above 
the accretion disk (with z(x+) > 0). An example for the mag- 
netic flux distribution is $(a;) = (x/x*)" 1 / 2 , resulting in mag- 
netic flux function <£>(.t) decreasing with radius faster than a 
monopole where = 1. 

Prescribing both the flux distribution and the shape of the 
flux surface does not over-determine the problem. The mag- 
netic flux function $ describes the opening of the magnetic 
flux tubes. With z(x), the shape of the flux surface chosen, 
the choice of the flux function just defines the position of 
the "other" flux surfaces. In a fully self-consistent approach, 
the field structure is determined by the solution of the Grad- 
Shafranov equation. Such solutions are not yet available. 



3.4.3. The plasma temperature 

The temperature distribution along the field line follows a poly- 
tropic gas law, T = (n/n*) r_1 . In our approach the tem- 
perature at the injection radius x* is determined by choosing 
the sound speed at this point, c s *, 



T* = 



r - 1 ( cL \ m p c 2 



(9) 



r \v-\-ciJ k B 

For typical parameters applied in our calculations, c s * = 0.05, 
r = 5/3 this gives a gas temperature of the disk corona of 
about 1.5 x 10 10 K at a jet injection radius x± = 8.3. This 
temperature is in rough agreement with the disk temperature of 
the advection dominated accretion disk models at small radii 
(Narayan et al. 1998). A smaller x+ requires a higher sound 
speed parameter implying a higher temperature T*. 

3.4.4. The iso-rotation parameter Op 

The iso-rotation parameter ^f^) of the field line is deter- 
mined from the position of the injection radius x*. This choice 
corresponds to the interpretation often applied for fip as the 
"angular rotation of the field lines". Here, we assume that the 
field lines are anchored in a Keplerian disk, f2p — ^Disk — 
17k (x*). The angular velocity of the last stable circular orbit 
around a Kerr black hole is Op(x*) ~ ±(x^ 2 ± a) -1 (the ± 
stands for co-rotation or retrograde rotation, respectively). For 
a radial position not too close to the black hole, the angular ve- 
locity in the accretion disk follows its Newtonian value. Close 
to a black hole ftp is limited due to the 'rotation of space' uj. 
An injection radius x+ = 8.3 gives £1f = 0.04 which is about 
0.1 for a = 0.8. 
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Fig. 3. Solution S3. Properties of the critical wind solution 
along a given flux surface (see parameters in Tab. 1). The small 
window shows the solution branches around the slow magne- 
tosonic point enlarged. The wind branch is the one with in- 
creasing velocity. The critical (magnetosonic) points are lo- 
cated at the intersections of the two solution branches (see Ap- 
pendix C for details). Top Poloidal velocity au p (in c). The 
asymptotic jet velocity of u p = 2.5 is reached after about 
x = 10 s . Middle Normalized proper particle density n (thick 
line) and temperature T in K (thin line). Below Normalized 
poloidal (thick line) and toroidal (thin lines) field strength, 
Bp, B^. Note that the injection radius is a;* =8.3. 

4. Results and discussion 

We now discuss our numerical solutions of the general rela- 
tivistic magnetic wind equation for different field geometries 
and input parameters. With the prescribed poloidal field our 
solution is uniquely defined by the conditions along the jet foot 
point and the condition of regularity across the magnetosonic 
points. Due to the stationarity assumption and the prescription 
of the field distribution, the spatial range of the computation is 
in principle not limited in radius. This is essential if one con- 
siders the huge size of Galactic superluminal jets in terms of 
the size of the central object. 

In general, we show that the acceleration of plasma from 
regions close to a black hole to the speed of 0.92 c observed for 
Galactic superluminal motion is possible to achieve. Depend- 
ing on the poloidal magnetic field distribution, the asymptotic 
speed of the jet is reached at a radius of about 100 gravitational 
radii. 

For comparison the leading parameters for our astrophys- 
ical solutions are summarized in Tab. 1. For illustration, we 
show the example solution SI demonstrating the typical fea- 
tures of the wind solution branches in the case of super- or sub- 
critical parameters (Fig. C. 1, Appendix C). The meaning of our 
figures is explained in detail in Appendix C. 

4.1. The wind solution - a collimating relativistic jet 

The time scale for the superluminal GRS 1915+105 jet is at 
least one month until the blobs become invisible in radio light. 
Mirabel & Rodriguez (1994) estimated that the ejection event 
for a blob lasts about 3 days. This time period would corre- 
spond to a value of Vt-p = 0.016 (for M = 5 Mq) and an 
injection radius of about x* ~ 15. The orbital period of the 
foot points rotating at the marginally stable orbit (for a = 0.8) 
is an order of magnitude less. The time scale derived for the 
intervals between the emission of jet knots is much larger as 
the period of the marginally stable orbit. The true location of 
the jet origin not yet known. Therefore, we suggest that the jet 
foot point should be located outside the marginally stable or- 
bit in order to maintain a jet flow for some time. For the first 
set of solutions we chose a foot point radius of = 8.3 or 
x+ = 15.3. 



The fact that the kinematic time scale of the blobs is at least 
10 times larger than the time scale for the generation of the 
blobs supports the assumption of stationarity in our calcula- 
tions. Clearly, on the long-term time evolution the presence of 
the blobs them self tells us that the jet flow is time-dependent. 

Compared to the other solutions in this sample with x* = 
8.3, solution S3 is weakly magnetized (Fig. 3). The initial open- 
ing angle of the magnetic flux surface is large (Fig. 2). The 
magnetic flux function $>(x) is constant along the field line. 
The asymptotic poloidal velocity of u p = 2.5 is reached be- 
yond a radius x ~ 10 8 (corresponding to a distance from the 
black hole of z{x) ~ 4 x 10 8 ). 

Figure 3 also shows the distribution of other dynamical 
variables. The poloidal field strength B p decreases with the 
opening of the magnetic flux surfaces. While the poloidal field 
distribution is prescribed in our approach, the toroidal magnetic 
field profile is a result of computation and therefore determined 
by the critical wind solution. At the injection point the toroidal 
field strength is about two times smaller than the poloidal com- 
ponent. Outside the Alfven radius the toroidal field becomes 
much larger than the poloidal component. For large radii the 
magnetic field helix is dominated by the toroidal component. 
In this region we find the toroidal field component following 
a power law distribution d(logS^) ~ d(\ogx). Therefore, in 
the asymptotic part the poloidal electric current is almost con- 
stant I ~ xB^ ~ const. In relativistic MHD electric — fields 
cannot be neglected. The electric field orientation is perpen- 
dicular to the magnetic flux surfaces and the field strength is 
\E±] = (R/ Rl)Bp. Therefore, the electric field is dominat- 
ing the poloidal magnetic field outside the light cylinder. 

Density and temperature are interrelated by the polytropic 
gas law. At the injection point the gas temperature T ~ 10 10 K 
(Fig. 3). The proper particle density at the injection point 
depends from the choice of the mass flux (in units of the mag- 
netic flux). Therefore, the calculated density profile n(x) may 
be applied to different mass flow rates (as long as the mag- 
netization a* is the same). Density and temperature decrease 
rapidly along the field line following the polytropic expansion. 
For x ^30 the proper particle density follows a power law 
rt/n* = 4 x 10 -5 a; -1 " 8 . At x ~ 1000 the gas temperature 
is about 10 6 K. Therefore we can estimate the size of a X-ray 
emitting region of about several 1000 r g in diameter. For the 
example of GRS 1915+105 this corresponds to 3.5 x 10 -9 arc- 
sec. It would be interesting to calculate the X-ray spectra of 
such an relativistically expanding high temperature gas distri- 
bution. 

Solution S3c2 has the same distribution of the magnetic 
flux function $ as solution S3. The magnetic flux surfaces, 
however, are collimating more rapidly. The derived critical 
wind solution has a higher magnetization, although the termi- 
nal speed and the total energy density E(^) of the S3c2 solu- 
tion is similar to S3. Because of the higher magnetization type 
S3c2 jet solutions have a correspondingly lower mass flow rate. 
The asymptotic speed is reached already at about x = 1000 
equivalent to a distance from the central black hole of about 
z = 3200. 
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Solution S9 relies on the same magnetic flux surface as 
S3c2. As a difference to S3c2, the magnetic flux function de- 
creases with radius implying a (spatially) faster magnetic field 
decay. As a consequence, the jet reaches its asymptotic velocity 
of Up = 2.57 even at about x = 100. The derived flow magne- 
tization is higher compared to S3c2 and S3 balancing the fast 
decay of the magnetic field distribution and we obtain the same 
asymptotic speed. This is interesting because it proves that not 
only the magnetization, but also the distribution of the magnetic 
flux along the field line determines the asymptotic speed. 

Note that solution S9 reaches the same asymptotic speed 
as S3c2 only because of its higher magnetization. Indeed, 
a solution similar to S9, but having the same magnetization 
er* = 1356 as for S3c2, only reaches an asymptotic speed 
of Up = 1.81 (not shown). Also, such a solution would be 
only very weakly magnetized in the asymptotic regime as the 
normalized flow magnetization changes as a ~ 1/ ^Jx for the 
$ ~ 1 solutions or a ~ 1/x for the $ ~ \\fx solutions 0, re- 
spectively. Similarly, in comparison, the asymptotical toroidal 
magnetic field is weaker by some orders of magnitude (a factor 
ten at x = 1000). In all the solutions presented in this paper the 
asymptotic jet is dominated by the kinetic energy. For the solu- 
tions with the large injection radius x+ = 8.3, the magnetic en- 
ergy is being converted into kinetic energy almost completely 
already at a radius of about several 100 gravitational radii. 

Solution S4 has the same magnetic flux distribution as 
S9, however, the field line is only weakly collimating. The 
asymptotic jet speed and the magnetization parameter is about 
the same. Only, the initial acceleration is weaker because the 
magneto-centrifugal mechanism works less efficient in the field 
with a smaller opening angle. 

Solution S4b has essentially the same field distribution 
as S4, but the injection radius is chosen larger. Therefore, 
the iso-rotation parameter fip is decreased by a factor of 
(8. 3/15. 3) 3 / 2 . As a result, a critical wind solution with a com- 
parable asymptotic speed could be obtained only for a very 
high plasma magnetization. This proves that highly relativis- 
tic jets can be expected even if the jet is not emerging from 
a region close to the black hole. Such a solution is feasible if 
the mass flow rate in the jet decreases with radius faster than 
the field strength (or flux distribution). The question remains 
whether such field strengths can be found at this position. 

We summarize the results of this section. The asymptotic 
speed is determined by the plasma magnetization and the dis- 
tribution of the magnetic flux along the field line. The shape of 
the magnetic flux surface determines the velocity profile along 
the field, thus, the position where the asymptotic velocity is 
reached. Highly relativistic outflows can be obtained even if 
the jet foot point is not very close to the black hole. However, 
in this case a high plasma magnetization is necessary. But this 

2 However, in the hot wind equation it is not possible to change 
only one single parameter in order to obtain a new set of critical wind 
solutions. In the case discussed above, with the decreased magneti- 
zation (i.e. an increased mass flow rate), the Alfven radius is corre- 
spondingly smaller (here, xa = 21.11 compared to xa = 22.93 for 
S9). 



seems to be in contradiction to the accretion disk theory (see 
below). 

4.2. The role of the magnetization 

The magnetic acceleration of jets and winds can be understood 
either as a consequence of converting Poynting flux (magnetic 
energy) to kinetic energy or due to Lorentz forces along the 
poloidal field line. In general, the higher the plasma magneti- 
zation the more energy can be transformed into kinetic energy 
of the wind. It has been shown theoretically for a cold wind 
that the relation between magnetization and asymptotic veloc- 
ity is that of a power law, itpoo ~ ctm 1 ^ 3 , for conical outflows 
(Michel 1969) and for collimating flows (Fendt & Camenzind 
1996), if the flux distribution is the same, respectively. How- 
ever, both papers do not consider gravity (and no general rela- 
tivistic effects). The new solutions presented in this paper are in 
general agreement with those results in the sense that a higher 
magnetization leads to a higher velocity. However, we are deal- 
ing with the hot wind equation and cannot derive a power law 
distribution from Tab. 1 , since the other wind parameters may 
vary between the different solutions. In difference to the cold 
wind solutions the magnetization is not a free parameter. In- 
stead, it is fixed by the regularity condition at the magnetosonic 
points. 

The wind magnetization is determined by the disk proper- 
ties at the jet injection points along the disk surface. For a stan- 
dard thin disk model that the ratio of the mass flow rate in the jet 
to the disk accretion rate is about 1% (Ferreira 1997). The ob- 
servational data for various jet-disk systems are consistent with 
this theoretical value. The accretion disk magnetic flux can be 
estimated assuming equipartition between magnetic field pres- 
sure (energy) and gas pressure (thermal energy) in the disk (see 
Sect. 3.4. 1). From Eq. (7) we find an equipartition field strength 
of about B cq ~ 5 x 10 8 G, if a V i S ~ 0.1 and R± = 10r g . Equa- 
tion (8) then defines an upper limit for the plasma magnetiza- 
tion at the injection radius, = 5 x 10 4 , for Mj et ~ 0.1M acc . 
Such a value is in general agreement with our solutions (Tab. 1). 
The maximum equipartition field strength estimated with the 
above given formulae can be much larger for Galactic black 
hole jet sources as for AGN (see Eq. 7). For a low black hole 
mass (with a smaller horizon) the disk comes closer to the sin- 
gularity and therefore becomes hotter. 

Again, we note that our estimate for the magnetization 
comes from comparison of different disk models (Sect. 3.4.1.). 
However, this does not mean that we apply a certain disk model 
for our computations. 

Finally, we come back to the wind solutions S4 and S4b. 
As already mentioned, these solutions demonstrate that the jet 
origin must not be necessarily close to the black hole. One may 
think that a strong magnetization at larger disk radii would do 
the job. On the other hand, the equipartition field strength in the 
disk decreases with radius implying that the highest magneti- 
zation and, thus, jet velocities must be expected from the inner 
part of the disk. Only, if the mass transfer rate from the disk 
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Fig. 4. Wind solutions S9 (upper left), S3c2 (lower left), S4 (upper right), S4b (lower right). Branches of poloidal velocity 
au p along the field line in units of the speed of light. For the solution parameters see Tab. 1. See caption of Fig. 3 for further 
explanation. 



into the jet decreases more rapidly with radius than the field 
strength, the magnetization increases. 

4.3. The influence of the rotating black hole 

As the main issue of our paper is the search for MHD wind 
solutions in Kerr metric, it is necessary to clarify the role of 
general relativity for the jet acceleration itself. Clearly, at an 
injection radius of i?* = 8.3 general relativistic effects are not 
very dominant. 

For comparison we have calculated wind solutions for a 
smaller injection radius x* = 3.3 (solution S3q, S3u, Fig. Dl, 
D2). The main effect is a much higher asymptotic velocity re- 
sulting from the rapid rotation, £l-p, at the smaller radius x+. 
With our choice Sip = 0.14 the asymptotic velocity drastically 
increases from u p = 2.503 (S3) to u p = 8.4792 (S3q). In or- 
der to obtain the critical solution for the higher rotation rate, 
the wind parameters have to be changed accordingly. o+ is de- 
creased by a factor of two, while c s * and it p * must be increased 
substantially. The large sound speed is in agreement with the 
smaller injection radius, since a higher disk temperature and 
pressure is expected close to the hole. The Alfven radius is de- 
creased by a factor of four, however, its location relative to the 
outer light cylinder remains the same. 

The limiting case of Minkowski metric can be achieved by 
setting M = and a = in the Boyer-Lindquist parame- 
ters (see Appendix A). For such a wind solution (S3u2) the 
magnetization is lower, although the asymptotic wind speed is 
the same as in the Schwarzschild case (see Fig. D2). This be- 
comes clear if we take into account that for S3u2 the wind flow 
does not have to overcome the gravitational potential Thus, less 
magnetic energy is needed to obtain the same asymptotic speed 
by magnetic acceleration. Further, we find from solutions for 
different angular momentum parameters a that in general the 
wind flow originating from a black hole with a smaller a is 
faster. As an extreme example we show the solution S3u3 cal- 
culated with a ~ but otherwise the same parameter set (see 
Fig. D2). This solution is magnetized stronger compared to the 
case of a = 0.8, thus, resulting in a higher asymptotic wind 
velocity. We believe that the reason for such a behavior is the 
fact that the effective potential of a black hole weakens (at this 
location) for increasing values of a. Therefore, less magnetic 
energy is necessary to overcome the effective potential. 

In the end, the results of this section are not surprising. 
They demonstrate that the wind/jet is basically magnetically 
driven. As a consequence, the acceleration takes place predom- 
inantly across the Alfven point as expected from MHD theory. 
Therefore, the scenario is similar to the case of classical pul- 
sar theory in Minkowski metric. For relativistic jets with a high 
magnetization the Alfven point is always very close to the light 



surface, which is defined by the angular velocity of the field 
line foot point. Usually, the Alfven point is located at a radius 
large compared to the gravitational radius. Thus, the influence 
of the general relativistic metric is marginal. Only, if the Alfven 
radius comes close to the hole, the choice of the metric will de- 
termine the jet acceleration. 

4.4. The question of collimation 

The huge size observed for the knots of the Galactic super- 
luminal sources leaves the possibility that the jet is basically 
un-collimated. 

Our numerical solutions have shown that the asymptotic 
speed of the jet does not depend very much on the degree of 
collimation in the flow. That speed is reached within a distance 
of about 10 8 r g . However, the observed upper limit for the knot 
size is still a factor 10 larger. Therefore, from our solutions, 
the observed knots are consistent with both a collimated and an 
un-collimated jet flow. In particular, solution S9 which is more 
collimated, has the same asymptotic speed as solution S4. 

In the case of extragalactic jets a high degree of collimation 
is indicated. The "lighthouse model" by Camenzind & Krock- 
enberger (1992) gives opening angles of only 0? 1 for the quasar 
2C273 or 0?05 for typical BLLac objects. The question arises 
whether there could be an intrinsic difference between the jets 
of AGN and Galactic high energy sources. Why should Galac- 
tic superluminal jets be un-collimated? A difference in the jet 
magnetization seems to be unlikely since the jet velocities are 
comparable. We hypothesize that if the jets of these sources 
are systematically different, this should rather be caused by 
the conditions in the jet environment. If the jets are collimated 
by external pressure, a different external/internal pressure ra- 
tio will affect the degree of jet collimation. Extragalactic jets 
are believed to be confined by an external medium (see Fabian 
& Rees 1995, Ferrari et al. 1996). It is likely that Galactic su- 
perluminal sources provide an example where the jet pressure 
exceeds the pressure of the ambient medium. While AGN jets 
bore a funnel through the galactic bulge, Galactic superluminal 
jets freely expand into the empty space. Such a picture seems 
to be supported by the fact that the Galactic superluminal jet 
knots move with constant velocity over a long distance. 

5. Summary 

We have investigated magnetically driven superluminal jets 
originating from a region close to a rotating black hole. The sta- 
tionary, general relativistic, magnetohydrodynamic wind equa- 
tion along collimating magnetic flux surfaces was solved nu- 
merically. The wind solutions were normalized to parameters 
typical for Galactic superluminal sources. 
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The assumption of stationarity allows us to calculate the jet 
velocity on a global scale over a huge radial range in terms of 
radius of the central source. The wind is launched close to the 
rotating black hole at several gravitational radii. The calcula- 
tion was performed up to a radius of 10 4 gravitational radii, but 
is in general not limited in radius. In some cases the asymp- 
totic speed may be reached only at a distance of several 10 8 
gravitational radii. Different magnetic field geometries were in- 
vestigated. The model allows for a choice of the shape of the 
magnetic flux surface and the flux distribution of that field. 

The physical wind solution is defined by the regularity con- 
dition at the magnetosonic points. As the poloidal field is pre- 
scribed, the choice of the following input parameters deter- 
mines the wind solution completely, (i) the injection radius of 
the matter into the jet, (ii) the injection velocity and (iii) the 
plasma magnetization (the ratio of magnetic flux to mass flux). 
The results of our numerical computation are the following. 

- In general, the observed speed for Galactic superluminal 
sources of more than 0.9 c can be achieved. 

- The flow acceleration is magnetohydrodynamic and takes 
place predominantly around the Alfven point. General rel- 
ativistic effects are important only if the wind originates 
very close to the black hole. In order to overcome the grav- 
itational potential, the critical wind solution must be higher 
magnetized in order to reach a similar asymptotic speed. 
This has been proven by calculating the Schwarzschild and 
Minkowski limit of the wind equation. 

- For a fixed magnetic field distribution the asymptotic jet 
velocity depends mainly on the plasma magnetization, in 
agreement with earlier papers (Michel 1969, Fendt & Ca- 
menzind 1996). The higher the plasma magnetization, the 
higher the final speed. The velocity distribution along the 
magnetic field shows a saturating profile depending on the 
distribution of the magnetic flux. 

- The magnetic flux distribution along the field line also in- 
fluences the plasma acceleration. Since the real field distri- 
bution is not known, we have considered two cases which 
show the typical behavior and which are probably close to 
the reality. We find that the jet velocity in a (spatially) faster 
decaying field can be the same as long as the magnetization 
at the injection point is high enough in order to balance the 
effect of the decrease in field strength. 

- For jet solutions not emerging from a region close to the 
black hole, a highly relativistic velocity can be obtained if 
the flow magnetization is sufficiently large. However, one 
we hypothesize that the field strength required for such 
a magnetization can be generated only close to the black 
hole. 

- Investigation of flux surfaces with a different degree of col- 
limation has shown that both field distributions allow for a 
relativistic velocity. The asymptotic jet velocity is reached 
considerably earlier in the case of the faster collimating flux 
surface. The jet reaches its asymptotic speed at a distance 
from the injection point of 3000r g or 10 5 r g , depending on 
the degree of collimation. The latter we measure with the 



opening angle of the collimating flux surface at this point 
and is about 15°or 45°, respectively. This distance is be- 
low the observational resolution by several orders of mag- 
nitude. Therefore, the question of the degree of collimation 
for the superluminal jets of GRS 1915+105 andGRO 1655- 
40 could not be answered. 

- Motivated by the huge size of the observed knots in the 
Galactic superluminal jets, we point out the possibility that 
the jet collimation process in these sources may be intrin- 
sically different in comparison to the AGN. For example, 
the upper limit for the knot diameter in GRS 1915+105 is 
about 10 9 Schwarzschild radii, which is distinct from typ- 
ical estimates for AGN jets with diameters of about 100 - 
1000 Schwarzschild radii. 

- The gas temperature at the injection point is about 10 10 K 
which is one order larger than the disk temperature at 
this point. With the polytropic expansion the temperature 
decreases rapidly to about 10 6 K at a distance of 5000 
Schwarzschild radii from the source. Both the temperature 
and the mass density follow a power law distribution with 
the radius. 

- The calculations show that the jet magnetic field is dom- 
inated by the toroidal component. Similarly, the velocity 
field is dominated by the poloidal component. 

In summary, our numerical calculations have shown that 
the highly relativistic speed observed for galactic superluminal 
sources can be achieved by magnetic acceleration. For a given 
magnetic flux surface we obtain the complete set of physical 
parameters for the jet flow. The calculated temperature, density 
and velocity profile along the jet would provide a interesting 
set of input parameters for computing the spectral energy dis- 
tribution. 
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for discussions and the analytical proof of the equations in Appendix 
B. We acknowledge fruitful comments by the referee Dr. Koide help- 
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Appendix A: Parameters of Kerr metric 

For the reason of completeness, here we list the parameters ap- 
plied in the equations of Kerr geometry. In Boyer-Lindquist 
coordinates with the parameters 

p 2 = r 2 + a 2 cos 2 9, A = r 2 + a 2 - 2M r, 

E 2 = (r 2 + a 2 ) 2 -a 2 Asin 2 6», Q = (E/p) sin 9, 

uj = 2aMrjcY?, a = pVA/H, 

the components of the metric tensor are defined as 

.goo = a m (2r/p(r,6) 2 - 1) 

.903 = -<r m 2rasin(9) 2 /p(r, 9) 2 

9n = a nl p{r,9) 2 /A(r,9) 

922 = a m p{r,9) 2 

.933 = <T ro E(r,0) 2 sto(0) a /p(r,0) a 

.9 = Det(g M „) = -.911522(530-500333) 
In our paper we have chosen a negative sign of the metric, 

Cm = -1. 
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Appendix B: Wind polynomial 

Here we provide the polynomial coefficients of the general rel- 
ativistic magnetohydrodynamic wind equation Eq. (4). For a 
derivation, see Camenzind (1986), Takahashi et al. (1990), or 
Jensen (1997). The specific angular momentum, properly nor- 
malized, is 

L = -(.903 + ^Fff33)/ (#00 + ^F.903) (B.l) 

For convenience we define the following parameters, 

r c °* („ P^^Y' 1 r r- n ^ 

Di = .900 + 2f2 F .9o 3 + ^f333, D 2 = (1 - n F L) 2 

D 3 = -(533 + 2Z.903 + L 2 g Q0 )/ (,9o3 - 300533) 

With the corresponding values at the injection radius x+ the 

total specific energy density of the flow E is defined as 

E 2 = -^« + l)(A, + ^ m M, 2 ) 2 

(D u + 2a m M 2 )D 2i , + D 3 *M$ ' 
where denotes the Alfven Mach number at the injection 
radius. The polynomial coefficients of the wind equation (4) 
are 



<Zl,2n+2m 




C l 2 


0-2,2n+m 




2a m C 2 Di 


^3, 2n 




D\ + C\ + (T m E 2 ClD 3 


0'4,2n—m 




2<r m C 2 D 1 + 2E 2 C 2 D 2 


a 5,2n-2m 




D\ + a m E 2 D 1 D 2 


^6,n+3m 




4C*iC 2 2 


Oi7,n+2m 






0>8,n+m 




2CiD\ + 4CiC| + a m 2E 2 dClD, 


0>9,n 




tr m 6C i .C 2 D l + 2E 2 dC 2 D 2 


fll0,n-m 




2C X D\ 


Oil, 4m 




6C 2 Cl 


Ol2,3m 




§o m C 2 C 2 D x 






C\D\ + 6CfCl + <j m E 2 C 2 ClD 3 


Ol4,m 




^C 2 C 2 D X 


«15,0 




C\D\ 


^16, 5m — n 




4CfC 2 


Ql7,4m-n 




2a IO C 1 C 2 D 1 


fll8,3m-n 






ai9,2m-n 




2o ra C\C 2 D\ 


O20,6m-2n 




etc 2 


a21,4m-2n 




etc 2 



All coefficients with the same second index have to be summed 
up, Ai = J2j - The polytropic indices n = 5, m = 3 give a 
polynomial of 16th order. 

Appendix C: Example wind solution in Kerr metric 

Here we show an example solution of the wind equation (4). 
The parameters are chosen such that a variation of <7* and c s * 
clearly demonstrates the criticality of the wind solution. They 
do not necessarily match astrophysical constraints. However, 



Fig. C.l. Example solution SI. Overlay of solutions u p (x) for 
three different parameter sets. 0+ = 49830, c s * = 0.4585 gives 
the critical solution which is regular across the magnetosonic 
points. The critical wind solution is the continuous branch start- 
ing with low velocity and accelerating to high speed. The mag- 
netization cr* is the critical parameter for the FM point, whereas 
c SVr is the critical parameter for the SM point. Sub- or super- 
critical solutions are obtained by variation of the parameters 
cr*, c s *. The choice of cr* = 51830, c s * = 0.4485 results in 
gaps in x(u p ), the choice of cr* = 48830, c s * = 0.4685 in gaps 
in u p (x). The other parameters are xa = 31.2, = 0.01, 
= 3.0, n F = 0.1 fi H = 0.025, a = 0.8. 

the asymptotic poloidal velocity is comparable to the speed of 
the Galactic superluminal sources. The solution (solution SI) 
considers a highly magnetized plasma flow with <r+ ~ 5 10 4 . 
The flux geometry is that of a slightly collimating cone with an 
opening angle decreasing with distance from the source. 

Figure CI shows the solution branches with a positive u 2 . 
An overlay of solutions for three parameter sets is displayed in 
order to show the typical behavior of wind solution. There is 
only one unique solution, the critical solution, with one branch 
continuing from small to large radii without any gaps in u p or 
x. The magnetosonic points are located at the intersections of 
the solution branches of the critical solution. The critical wind 
solution is regular at all three magnetosonic points. It is de- 
fined by a unique set of the parameters E, L and a (for flp pre- 
scribed). In the critical solution the slow magnetosonic point is 
passed close to the foot point of the jet. The Alfven point is lo- 
cated at x — 31 and the fast magnetosonic point not far beyond. 
The asymptotic speed of the flow is u p = 2.28, equivalent to 
v p ~= 0.9 c (not shown in the Figure). 

Sub- or super-critical wind parameters lead to solution 
branches which are not defined for all radii or all velocities. 
Even for a slight variation of these parameters the solution 
will be not continuous anymore, implying "jumps" or "shocks" 
across the gaps in the solution branches. At these locations the 
stationary character of the solution most probably breaks down. 
Such solution branches are inconsistent with the assumptions 
and are therefore referred to as unphysical. 

Appendix D: The wind solution for a small injection radius 

For comparison, we show solutions of the wind equation with a 
small injection radius X-^. — - 3.3 3.s well as solutions in the limit 
of Minkowski and Schwarzschild metric (for a discussion see 
Sect. 4.3). Solution S3q corresponds to solution S3, however, 
with a magnetization smaller by a factor of two. The asymptotic 
speed is u p = 8.48 and much larger than for S3. Also solutions 
S3u, S3u2, S3u3 correspond to S3 and S3q. However, in this 
case the Alfven radius and the derived magnetization parameter 
are lower resulting in a lower asymptotic speed. Solution S3u is 
the Kerr solution for a = 0.8, S3u3 the Schwarzschild solution 
(a = 10~ 8 ), and S3u2 the Minkowski solution where we set 
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Fig. D.l. Example solutions with a small injection radius X* = 
3.3. Rotation rate Hp = 0.14. Solution S3q with a = 0.8, xa = 
5.83, o~ * = 480 has a high asymptotic velocity u p ~ 8.48. 

Fig. D.2. Example solutions with a small injection radius x* = 
3.3. Rotation rate fi F = 0.14. Top Solution S3u with a = 0.8, 
xa = 5.33, and a smaller <7* = 100 and a lower asymptotic 
velocity u p = 2.96. Middle Solution S3u3 in Schwarzschild 
metric, a = 10~ 8 , xa = 5. 33. The asymptotic speed is u p = 
4.48 with cr* = 205.7. Bottom Solution S3u2 in Minkowski 
metric, a = 10~ 8 , M = 0. The asymptotic speed is u p = 4.55 
while er* = 82.5. 

a = 10~ 8 and M = in the Boyer-Lindquist parameters (see 
Appendix A). For a comparison of all solutions see Tab. 1. 
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